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We study the use of entanglement purification for quantum communication over long distances. 
For distances much longer than the coherence length of a corresponding noisy quantum channel, 
the fidelity of transmission is usually so low that Standard purification methods are not applicable. 
It is however possible to divide the channel into shorter segments that are purified separately and 
then connected by the method of entanglement swapping. This method can be much more efficient 
than schemes based on quantum error correction, as it makes explicit use of two-way classical 
communication. An important question is how the noise, introduced by imperfect local operations 
OS ' (that constitute the protocols of purification and the entanglement swapping), accumulates in such 

a compound channel, and how it can be kept below a certain noise level. To treat this problem, we 
first study the applicability and the efficiency of entanglement purification protocols in the situation 
bX), of imperfect local operations. We then present a scheme that allows entanglement purification over 

arbitrary long channels and tolerates errors on the per-cent level. It requires a polynomial overhead 
' in time, and an overhead in local resources that grows only logarithmically with the length of the 

channel. 
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In quantum communication, qubits jl| are sent across quantum channels that connect distant nodes of a quantum 
network. In general, these quantum channels are noisy and therefore limit the fidelity of the transmission. A bottleneck 
for quantum communication is the scaling of the probability of transmission errors with the length of the channels 
connecting the nodes. Different from classical (digital) communication, the "signals" in quantum communication 
consist of single qubits which may be entangled. Owing to fundamental principies, these qubits cannot be cloned ^|J|] 
Qh, nor amplified Q without destroying the essential quantum feature. This fact limits the maximum distance between 
+-> ' the nodes to a few múltiples of the absorption length (or the coherence length) of the channel, and poses a severe 
restriction on any practical application. For example, in recent experiments on quantum cryptography single 
photons are sent through optical fibers with a given absorption and depolarization length. This has two effects: (i) 
to transmit a photon without absorption, the number of trials scales exponentially with l; (ii) even when a photon 
• • i arrives, the fidelity of the transmitted state decreases exponentially with l. In the experiment, the distance between 
. ^ ' the nodes is therefore presently limited by a few múltiples of the absorption length of the fiber . 

In the context of fault-tolerant quantum computing |?J , Knill and Laflamme |8j have discussed an important scheme 
that allows, in principle, to transmit a qubit over arbitrarily long distances with a polynomial overhead in the resources. 
Their method requires to encode a single qubit into a concatenated quantum code (i.e. an entangled state of a large 
number of qubits) and to operate on this code repeatedly while it is propagated through the channel. The tolerable 
error probabilities for transmission are less than 10 -2 , whereas for local operations they are less than 5 x 10 -5 . This 
seems to be outside the range of any practical implementation in the near future. A crucial figure for any experiment 
will be the number of partides that can be manipulated locally in a coherent fashion, together with the precision 
with which such local manipulations can be realized. For instance, even when the number of partides scales only 
polynomially with the length of the channel, a polynomially large number may still be far too large for any practical 
implementation. 

A central tool in the theory of quantum information, specifically for quantum communication is entanglement 
purification p|-|Ï2|], It allows, in principle, to create maximally entangled states of partides at different locations 
even if the channel that connects those locations is noisy fi]] . These entangled partides can then be used for faithful 
teleportation [^3| or secure quantum cryptography JÏ2|,|Ï4j. It seems therefore natural to use this method as an 
ingredient for a quantum repeater. Furthcrmore, it allows highly efficient two-way protocols which can not be realized 
with quantum error correction procedures. If a noisy quantum channel is much longer than its coherence length, one 
cannot directly employ Standard purification schemes. These schemes first create an ensemble of low fidelity EPR 
pairs across the channel, and then purify/distill a few perfect EPR pairs out of the ensemble. This, however, requires 
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a minimum fidelity F m i n of the initial pairs to operate with, which cannot be achieved as the length l of the channel 
increases. 

In this paper, we treat the problem of a quantum repeater based on the method of entanglement purification. The 
idea of such a repeater is to divide a long quantum channel into shorter segments, which are purified separately, before 
they are connectcd. Connecting two segments of a channel means here to build up quantum correlations across the 
compound channel from correlations that exist across the individual segments. This can be done by entanglement 
swapping, or teleportation of entanglement. A quantum repeater must therefore combine the methods of entanglement 
purfication and teleportation. 

Although the combination of these methods should, in principle, allow to create entanglement over arbitrary dis- 
tances, it is another question how much this "costs" in terms of resources needed for purification. Resources means 
here the number of low fidelity EPR pairs that have to be provided for purification of each channel segment. This 
quantity is related to the number of partides that have to be manipulated locally (at the connection points between 
the segments) in a coherent fashion. If the resources grow too fast with the length of the channel, not much will 
be gained by the whole procedure. A further important quantity is the error tolerance for the local operations. In 
every real situation, the local operations applicd to one or more partides, will bear some imperfections. Since such 
operations are the building blocks for any entanglement purification protocol, their imperfections will limit the max- 
imum attainable fidelity for an EPR pair and the efhciency of the protocol. In the context of the quantum repeater, 
a maximum fidelity F < 1 corresponds to a residual amount of noise for each segment. When the segments are 
connected, this noise aceumulates. 

To give a realistic treatment of the quantum repeater, we will therefore first study the applicability and the efficiency 
of entanglement purification protocols in the situation of imperfect local operations. The general conditions under 
which Standard purification protocols can be used in the presence of errors have been studied by Giedke et al. [ p"5[ . 
This includes, in particular, thresholds and lower bounds for the attainable fidelities. For a genèric class of stochastic 
errors, one can give explicit representations of the imperfect operations in terms of completely positive maps and/or 
POVMs JTsJ . In terms of these maps, we derive recursion formulas for the fidelities, which generalize the results given 
by Bennett et al. || and by Deutsch et al. ||ï^| . From these results, we estimate that, for a genèric class of errors, 
Standard entanglement purifeation protocols work even for error probabilities of the order of a few per cent. In the 
context of long distance communication, we develop a purification procedure that combines the Standard protocols with 
quantum teleportation in a certain way to be specified. This purification procedure allows quantum communication 
via noisy channels of arbitrary length. Since it explicitly exploits two- way classical communication 10 1, our procedure 
is much more efficient for quantum communication than protocols based on quantum error correction || . Specifically, 
this solution of the quantum repeater tolerates errors on the per-cent level; it requires a polynomial overhead in time 
and an overhead in local resources that grows only logarithmically with the length of the channel. 

The paper is organized as follows: In Sec. |ïï| we introduce a genèric error model to describe imperfect local operations 
such as noisy 1 bit and 2 bit operations and imperfect measurements. With the aid of these models, we revisit in 
Sec. III the bàsic elements of quantum communication - purfication and teleportation. This discussion includes results 
for the maximum attainable fidelities and the working conditions of the Standard protocols under the condition of 
noisy local operations. After this discussion, in Sec. IV, we are ready to attack the quantum repeater problem itself. 
Our solution combines the methods of entanglement purification and entanglement swapping into a single (meta-) 
protocol, which we call nested entanglement purification. This protocol allows to distribute entanglement with a 
given fidelity across arbitrary long distances, even when local operations are imperfect. It requires an overhead in 
local resources that grows polynomially with the length of the channel. A modification of this protocol works with 
resources that grow only logarithmically with the length of the channel (while the build- up time grows polynomially) . 



II. IMPERFECT LOCAL OPERATIONS AND MEASUREMENTS 



Imperfections of local operations can have various origins. An example for an imperfect one-qubit operation is a 
unitary 'overrotation' [|Ï7|], as it would be realized by a làser pulse that effects a Bloch rotation around an angle that 
slightly deviates from the ideal (intended) one. If the làser is well stabilized but the pulse area tuned incorrectly, 
say, the deviation of the real from the intended (unitary) rotation will be the same when we repeat the operation, 
representing an example of a systematic error. If, on the other hand, the angle of the rotations fluctuates randomly 
around a certain mean value, every pulse introduces a certain amount of noise into the system, and the resulting 
operation would be described by a non-unitary map. In the present paper, we concentrate on this latter situation, 
that is we assume only stochastic errors and neglect systematic errors. 

To give a precise description of a noisy 1-bit or 2-bit operation, one needs to know the exact form of the error 
mechanisms, which in turn depend on the specific physical implcmentation. In general, we have only limited knowledge 
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of these details. A genèric model for a noisy channel is the so-called (symmetric) depolarizing channel. It transforms 
a qubit with initial state p in the way 



P^PP+^I, (1) 

where / is the identity operator and < p < 1 depends on the time for which the qubit is subjected to the channel. 
The action of the channel thus results in admixing a completely depolarized state 1/2 to the initial density operator. 
The limit p — > corrcsponds to a very noisy channel, while p — > 1 describes a channel with very little noise ("ideal 
storage" ) . 

One can easily generalize ^ to the case where the ideal operation is some arbitrary unitary (1-bit) operation U, 

p^pUpUÏ + ^I. (2) 

In this situation, the ideal operation is accompanied by the action of the depolarizing channel. The whole map 
becomes non-unitary. It describes an imperfect 1-bit operation, which for p — > 1 becomes identical to the ideal 
unitary map. We call p the reliability of the imperfect operation, which represents a lower bound to the probability 
that its result corresponds to the ideal one. For a parametrization U — Us = exp(—ia ■ a) an imperfect (noisy) 
rotation may explicitly be defined by the map 



P — R s p = jdp v0)U s+0P UÍ +0 . (3) 

where the noise function is normalized to unity and describes the fluctuations Sa = (3 of the (generalized) Bloch 
vector a. One can easily show that for "isotropic noise," v(j3) = z>( |/3 1 ) , the integral reduces to 

RsP = pU s pUÍ + ^-^I 

= pR^p+i-P I (4) 



wherein R l ~p = Usplft describes the ideal rotation. 

Motivated by these observations, we generalize this model to arbitrary one-qubit and two-qubit gates |D| acting 
on an entangled state p of several qubits. We thereby require that errors introduced by local operations may only 
affect the local qubits on which the ideal operation acts (locality condition). Let p be a state of q qubits, labeled by 
1,2, ... ,q. An imperfect one-qubit operation acting on the first qubit, say, is then described by the map 

O lP = p 1 Of^p+ l —^tv l {p}®h (5) 

whereas an imperfect two-qubit operation acting on qubits 1 and 2 is described by 

12P = P2 Ofr 1 P + ^r^tn 2 {p} ® 1x2 . (6) 



In these expressions, O ldcal is the ideal (perfect) operation and 1\ and 1\ 2 denote unit operators on the subspace 
where the ideal operation acts, corresponding to totally depolarized one and two-qubit states. The state of the 
other partides are described by the partial traces tri{p} and tr\ 2 {p}, respectively, of the initial density operator 
over theses subspaces. The quantities p\ and p 2 measure the reliability of the operations, where perfect operations 
correspond to pj — 1 . Technically speaking, an imperfect operation is modeled by the corresponding ideal operation 
accompanied/followed by a depolarizing channel that acts (only) on the same subspace as the ideal operation. 

Note that the maps (gj) and (|6|) are linear and trace conserving, albeit non-unitary. For any state p which is diagonal 
in the Bell basis defined by any two of the partides - all states we deal with in this paper are of this form -, the model 
is self consistent in the sense that 0\ ® 2 (two single-qubit operations, each described by an error parameter pi) can 
be written as a certain two-qubit operation (9 12 described by an error parameter p 2 = p\. Generally, any sequence 
of local operations can be written as a single joint operation within this model. The resulting error parameter (i.e. 
reliability) describing the joint operation is obtained by multiplying the error parameters of the single operations. 

We will usc these maps to estimate the role of imperfect operations in Standard protocols for entanglement purifi- 
cation and for teleportation. Although the model is simple, it is genèric at the same time. We just mention here that 
similar maps may be derived from a Jaynes-type principle, insofar as their resulting states maximize the entropy of 
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the resulting state with respect to variations that are trace conserving and keep the average fidelity of the operation 
constant. 

We finally describe an ímperfect measurement on a single qubit (in the computational basis) by a POVM pc| ] 
corresponding to the 'imperfect projectors' 

Pg = V \0){0\ + (l-r))\l)(M 

iï=»,|l>(l| + (l-» J )|0)(0|. (7) 

The parameter r\ is a measure for the quality of the projection onto the basis states. Assume the qubit to be measured 
is in the state p = |0) (0 and we are trying to measure its state with the aid of an measurement apparatus described by 
(0). The expectation vàlues (Pq) = r\ and {P?) = 1 — r\ simply mean that the apparatus will give us the wrong result 
("1") with probability 1 — r\ > 0. That is, the result is not completely reliable. An ideal measurement, in contrast, is 
described by r\ — 1. It is clear that the effect of the measurement on the measured qubit is not fully specified by the 
POVM (Q). In the present context, however, this description is sufficient as all measured partides are removed from 
the system i.e. we trace over their degrees of freedom after the measurement. 

The operations (JsJ) and (||) together with single-qubit measurements (0) are sufficient to describe all operations 
occurring in the context of teleportation and entanglement purification. For example, a Bell measurement (the 
measurement of a projector in the Bell basis) on two partides, say 1 and 2, can be realized by a two-qubit operation 
0\2 = CNOT 1 ]™^ 01 (controlled NOT) followed by a Hadamard rotation of particle 1 and two single-qubit measurements 
Pç,(iy ^o(i) on P ar tides 1 and 2 |p5| . Instead of performing the Hadamard transformation one can also measure particle 
1 in a rotated basis. In summary, an imperfect Bell measurement is described by an imperfect two-qubit operation 
followed by 2 imperfect single-qubit measurements, effecting e.g. |0)|0) ± |1)|1) — > (|0) ± |1))|0). 



III. PURIFICATION AND TELEPORTATION WITH IMPERFECT MEANS 



In this section we reconsider the bàsic elements of quantum communication, which are teleportation and purification, 
in the presence of local errors. 



A. Entanglement purification 



Purification is the distillation of few 'perfect' EPR pairs out of many imperfect pairs. In the following, we will 
generalize two different purification protocols that have been treated in the literature, by introducing imperfect 
gate and measurement operations. In subsection III A 3 , we then discuss a modified version of these schemes which, 



generally speaking, has less favorable convergence properties but works with a smaller and constant number of physical 
resources. 



1. Scheme of Bennett et al. (Scheme A) 

This purification scheme was introduced by Bennett et al. ]9|,[Ï0|]. In short, the scheme takes two pairs (1-2 and 3-4) 
as in Fig. [|, both in a Werner state 

P = ^i$+)(<d+i + (^) (hn<in + + \*-)(*-\) ■ (8) 

with fidelity F — (<fi + \p\(f) + ) . Then it performs local (1 & 2-bit) operations on the partides at the same ends of the 
pairs. Writing pi2 ® p^i = P1234, this essentially involves two CNOT operations, CNOT' 1 ™ pcrf and CNOT™ 1301 ^ acting 
on the state P1234, followed by a simple measurement P^iy ^0(1) °^ ^ ne P ar ticles 3 and 4. If the partides are found in 
the same state (00 or 11), the remaining pair, described by the state p' 12 is kcpt, otherwise it is discarded. To obtain 
a recursive mapping between Werner states, a non-unitary depolarization operation (twirl) is applied to the resulting 
state before it is used for the next purification step. The original treatment of this protocol [p|JÏ0|] assumes that all 
operations and measurements of the protocol are perfect. 

To estimate the effect of local errors, we evaluate the protocol with the aid of the imperfect operations (^-0). The 
calculation of the fidelity F' of the new pair p' 12 is somewhat lengthy but straightforward, and the result can be given 
in analytic form. The result of this calculation is summarized in the formula 
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[F 2 + (^)W + (1 - V?] + [F(±f) + (^) 2 )} [2^(1 - „)] 



-P2 • 



F' = : _ L __^ : ■____ L _ *_ L _ : _ 2 - ggi ; (9) 

[F* + - F) + |(1 - F) 2 ][r] 2 + (1 - + [F(±f) + (^) 2 )]N(1 - 7?)] + 4(^) 

which reduces to the formula given in Ref. (|] in thc limiting case 77 = 1 and P2 = 1 (perfect operations). Although 
p' 12 is no longer a Werner state, it can again be brought to Werner form using depolarization, which is considered to 
be noiseless here ]Ï9fl . If one starts with an ensemble of pairs with fidelity F, F' gives the fidelity of the remaining 
pairs that are left after one purification step. This defines a fraction of 2/p cvcn of the initial ensemble of pairs, where 
Peven denotes the probability for finding the partides 3 and 4 in the state 00 or 11, and is given by the denominator 
in (JoJ) . The factor 2 arises since the pair originally composed of partides 3 and 4 is lost due to the measurement. 

Each (iterated) application of (||) corresponds to a purification step. For each purification step, 2 identically pairs 
(which both result from previous successful purification steps) are used. The purification resources are dcfined by the 
average number M of pairs needed to perform k max successful purification steps, and are given by 

£ Peven 

(k) 

where the probabilities Povcn depend on thc purification step k. The physical resources, defined by the total number 
of partides at A or B that are used to store the pairs needed for the purification process, are for this scheme equal to 
the purification resources, since all pairs are created at the beginning and have to be stored as can be seen in Fig.^. 

When calculating the fixpoints of (^|), one finds that one fixpoint is always F — j, independent of the error 
parameters pi and 77. The other two fixpoints -F m i n and F max (see Fig. |ï^) are given by the expression 



_ 877(77 - 1) + 3 ± a/10 - 9/pj + 64t7 4 - 128?? 3 + II677 2 - 52r/ - 3677(77 - \)jp\ ^ . 

~ 1677(77- 1) +4 ' ' ' 



They depend on the error parameters and give the borders of the interval within which purification is possible. If 
F 6 (F m i n , F max ), then F' > F. Since F max is an attractor, iterative application of (^) leads to a resulting pair with 
fidelity F — ► F max . The value F m i n thus gives the threshold for F where this purification scheme can be successfully 
applied, while F max gives the maximal reachable fidelity. For perfect local operations, F m [ n = | and F max = 1, 
meaning that all pairs with F > \ can be purified to F = F msüC = 1. For imperfect local operations, F m j n > | 
and F max < 1, i.e. no perfect EPR pairs can be created. If the error parameters become larger, these two fixpoints 
approach each other and the region where purification is possible shrinks to zero. The limiting situation f roax — -Fmin 
defines the threshold for the applicability of the purification protocol. For all pairs (j>2,v) f° r which there is only one 
real fixpoint (at F = 1/4), the imperfections of the local operations introduce more noise than one gains from the 
purification, so the scheme breaks down. For example for 77 = 1, (perfect measurements) , Eq. ( |ïl"| ) simplifies to 
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F °" =i ± iV 10 ~ 9 / p 2 ( 12 ) 



with the threshold at P2 = y/9/10 ~ 0.95. That is, the CNOT gate must work with a reliability of at least 95%. The 
threshold gets tighter when the measurements are imperfect as well, i.e. for 77 < 1. In general, the threshold vàlues 
for the error parameters P2 and 77 are of the order of some percent as can be seen in Fig.|3j. The fixpoints as a function 
of the error parameters are plotted in FigM where we set P2 — 77 to get a 2 dimensional plot. 



2. Scheme of Deutsch et al. (Scheme B) 

A purification protocol that converges faster and involves less resources was proposed by Deutsch et al. [ jÏ2| . 
Generally speaking, this scheme can be described as a mapping on states that are diagonal in the Bell basis, but need 
not necessarily be Werner states, 

pia = A\4>+){4>+\+ B\i>-){,p-\ + C\i> + )(i> + \ + D |0-><0-| (13) 

where |^ ± ) and ïp ± ) are the Bell states between partides 1 and 2 in the usual notation. In the notation of (|Ï3|), the 
protocol corresponds to a mapping R 4 — » i? 4 between the diagonal clements (A,B,C,D). In particular, between two 
successive purification steps, the states are not depolarized. This feature, together with additional n/2 rotations that 
are applied to the qubits before the CNOT operations, makes scheme B a much faster converging protocol. 
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In the following analysis, we evaluate scheme B using the imperfect operations (||)-(0) instead of perfect ones. A n/2 
rotation followed by a CNOT operation is thereby treated as a joint 2-qubit operation with a single error parameter 
P2- The resulting map between the diagonal elements i? 4 — > i? 4 is described in p2[ ; it completely characterizes the 
action of this purification protocol. Successive purification steps are now described by iterated applications of this 
map. For each purification step, two identical pairs (which result from the previous step) are used. The resources can 
be calculated in a similar way as for scheme A. 

The fixpoints of this map are no longer described by a single parameter F but by a set of 4 numbers 
(Afi X , Bfi X ,Cfix, D fi X ). The fidelity of such a diagonal Bell state is given by A, the \4> + ) component and, for 
simplicity, we shall continue to call this component "fixpoint" of scheme B, although the fixpoint is not sufficiently 
described by this single parameter. For example, it may happen that a Werner state with a certain fidelity Fq cannot 
be purified, while a binary state with the same fidelity F that has only two nonzero components \(f> + ) and say, 
(in short: "binary state (A,C)") can be purified. 

We have numerically compared schcmes A and B, and found that scheme B converges much faster towards the upper 
fixpoint F max if Fq > F min as can be seen in Fig.|| Surprisingly, the upper fixpoint for scheme B is above the one for 
scheme A for given error parameters and does not depend on the "shape" of the initial state. The lower fixpoint F m í n 
is smaller for scheme B and therefore the interval within which purification is possible is larger for scheme B, and so 
is the maximal reachable fidelity. The fixpoints as a function of the error parameters P2 = V are plotted in Fig.|^. One 
can see that scheme B is significantly less sensitive to noisy local operations than scheme A. 



3. A modified purification scheme (scheme C) 

Both schemes (A and B) need many elementary pairs to start with, which have to be stored somewhere before the 
protocols are applied. The physical resources are thus quite large and grow with the number of necessary purification 
steps. In the following, we modify scheme B in a certain way, and call this scheme C. 

Different from scheme B, we do not use two identical pairs at each purification step (which are left over from 
previous purification steps) but always purify one and the same pair with the help of an auxiliary pair 7r with 
constant fidelity A Vo . Apart from that, we employ the same protocol (that is the sequence of local operations and 
measurements) as in scheme B. The resulting map is thus the same, the only difference being that P34 is a diagonal 
Bell state which is constant throughout the whole purification procedure and is given by (A wo , B^ a , C^ , D no ) . This 
auxiliary pair is repeatedly created before each purification step between two of the four partides as visualized in 
Fig|| If a purification process is not successful, one has to start from the beginning (when both pairs 12 and 34 are 
put in the state ttq). Using this procedure, the fidelity converges towards a given fixpoint F(tto) which is not equal to 
one E] . The fixpoint highly depends on both the initial fidelity and the shape (other diagonal elements in the Bell 
basis) of the state ttq . The properties of the fixpoint will be discussed in more detail later on. 

For ideal operations, this method is generally less favorable as far as its convergence properties are concerned. 
However, for imperfect operations, the situation changes and the drawback that the reachable fixpoint F %0 is smaller 
than unity becomes less important. As we shall see later, for imperfect local operations the fixpoint of scheme C may 
lie even above the fixpoints of schemes A and B (see Fig. |J). The main advantage of this scheme is, however, that 
the physical resources, that is the number of partides needed at A or B to store the pairs, are constant (namely 2) 
and independent of the number of necessary purification steps. This fact makes scheme C particulary interesting for 
the repeater problem, where the aceumulation of local partides at the connection points of a compound long channel 
plays an important role. 

The purification resources, on the other hand, that is the (average) number of how often an auxiliary pair has to be 
created in order to perform k max purification steps, have now to be calculated in a different way. Let p e ven{k) be the 
probability to succeed at the k th purification step. The total resources can be calculated using the iteration formula 

M "+' = (M " + 1) ï3rrT)- (14 > 

To close this formula, note that Mq=1. This can be understood as follows: Starting at the elementary level, one needs 
two pairs for purification and has to repeat the procedure on average - — ^-^y times. Thus M\ = (1 + 1)- — ^jj • For 
the next purification step, one needs one additional elementary pair to perform the next purification step. This step 
has to be performed on average - — 3— times. The resources to perform the second purification step are thus given 

by M2 = (Mi + 1) — — í-^y and so on. To perform k max purification steps, one needs to create on average 
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times an auxiliary pair. 

A related quantity is the average number of purification steps (including also those which fail). This is given by 
the iteration formula 



S k+ i = (S k + l) (16) 

Peven < L ) 

where So = and k denotes the number of successful purification steps. If no purification is performed, the number of 
steps S'ò = 0. This is the only difference to the iteration formula (|Ï5| ) for the resources, since every time an elementary 
pair is created, one performs a purification step. Thus the average total number of purification steps performed 
(assuming that k max steps were successful) is given by 

g(n(^)). 

The fixpoint discussion for this purification protocol is more involved due to the many parameters, namely the 
fidelity A Vo and the shape (other diagonal elements B Vo , C wo , D^ )) of the elementary pair, and the error parameters 
P2 and rj describing the noisy local operations. For this protocol, purification is possible if the rcachable fidelity 
(fixpoint) for a certain auxiliary pair lies above its initial fidelity. To illustrate first the shape dependence of the 
fixpoint, we introduce the following parameterization of diagonal Bell states with constant fidelity Fq. In the notation 
of ([H^) we write 

A = F 
S = C=(l-F )(l-e)/2 

D = (l-F )e (18) 

where < e < 1. A Werner state thereby described by e = -| and a binary state (A,D) is obtained for e = 1. Figure 
[j] shows the dependence of the fixpoint as a function of e for different error parameters. 

For the manifold of states covered by this parametrization one sees that the 'optimum shape' of the elementary pair 
is a binary state (A,D), meaning that the reachable fixpoint for a given initial fidelity is the largest for this shape. 
Figure H demonstrates two different things. First, it shows that the fixpoint (maximally reachable fidelity) depends 
on the initial fidelity of the elementary pair even if the other components are kept constant. Second, it demonstrates 
that the reachable fidelity using this purification scheme may lie above the reachable fidelity using schemes A or B. 
(The upper fixpoint for scheme A is not plotted here, but as we already know (see Fig.|^) it always lies below the 
fixpoint for scheme B). One should however mention that this is only the case if the shape of the elementary pair is 
close to the optimal configuration (A,D). This can be understood as follows: Scheme B always picks 2 pairs that result 
from a previous purification step and the shape of such a pair converges towards a "working state" which no longer 
depends on the initial shape but only on the error parameters. Scheme C, on the other hand, always uses the same 
auxiliary pair for each purification step which has not been influenced by noisy local operations in previous steps. As 
we have pointed out earlier, in general one can only increase the fidelity by a certain amount. 



B. Teleportation 

In the context of the quantum repeater, teleportation comes into play when two (purified) segments of a channel 
are connected. The channel segments are represented by EPR pairs between partides 1-2 and 3-4 that have been 
created across the corresponding segments, see Fig. [| In general, to teleport the state of particle 2 to particle 4, a 
Bell measurement has to be made on partides 2 and 3 (and the result of this measurement communicated particle 
4). If particle 2 is itself entangled to some other particle, say particle 1 as in Fig. |^, the effect of the teleportation 
is to transfer this entanglement to an entanglement between partides 1 and 4. This process has also been termed 
"entanglement swapping" p4j . 

As a result, one obtains an EPR pair between partides 1 and 4 which can be used as a single channel (again with 
the help of teleportation). In this sense, the segments of the compound channel have been connected. Put in different 
terms, the connection process creates quantum correlations across the compound channel from correlations that exist 
across the individual segments. 

We will study this process first for the particular case where the involved pairs (1-2 and 3-4) are not maximally 
entangled but are both in a Werner state (||) with fidelity F. To connect the pairs (1-2) and (3-4), we make a Bell 
measurement on partides 2 and 3, which is realized by a (impcrfcct) CNOT operation followed by two (imperfect) 
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single-particle measurements, as described in the last paragraph of Section H. Depending on the outcome of this 
measurement, a one-particle operation O4 is performed on particle 4, as in the teleportation scheme |Ï3f| , The 
measurement outcome needs thereby to be sent to particle 4 using classical communication. The result of this 
sequence of operations, followed by a subsequent depolarization, is a Werner state between particle 1 and 4 with a 
smallcr fidelity F2, which is independent of the outcome of the measurement. 



More generally, imagine that we have not only 2 pairs but a whole string of N pairs as visualized in Fig. [10[ Each 
of the pairs (A — Ci),(Ci — C2),....(CW-i — B) is assumed to be in a Werner state with fidelity F. Connecting these 
N pairs using the procedure described above leads to a pair (A — B) with fidelity -F/v given by the formula 




3 / v 

The parameters rj,pi and P2 appearing in this formula quantify the amount of noise that is introduced by the 
connection processes. The connection therefore leads to an exponential decrease of the resulting fidelity, unless both 
the elementary pairs and all the operations involved have unity fidelity. 

The connection can be performed in two different ways, which both lead to the same resulting state with fidelity 
-Fjv but involve a different temporal ordering of the operations. The first way is to connect the pairs sequentially, i.e. 
first connect at Ci, then at C2, and so on, each time only connecting one additional pair. This involves a sequence 
of N — 1 connection procedures. The second way is to connect the pairs in parallel. To achieve this, first connect 
simultaneously the neighboring pairs at C\, C3, ... Cjv-i- This lcaves us with longer pairs (A — C2), (C2 — C4),..., 
(Cjv_2 — B). Then connect simultaneously these longer pairs at C2,Cq,...,Cn-2, and so on, until we get a final pair 
between A and B. To have at each step pairs of equal length and fidelity, N should be some power of 2, N = 2™, 
although this is not an essential requirement. This method is much faster as it requires fewer iteration steps, namely 
log2N = n instead of N — 1, although the number of local connections to be made is, of course, the same. 

The elementary pairs that are connected need not be in a Werner state. For example, for Bell diagonal states one 
can derive a map R A — > R A that relates the elements (A, B,C, D) of the corresponding states before and after the 
connection. That is, Bell diagonal states are mapped to Bell diagonal states under connection. All other remarks 
about the connection procedure for a string of pairs apply here as well. The parallel way of connecting the pairs 
is preferable not only because can be done faster, but also because, after each connection step, one deals with an 
ensemble of identical pairs. Further details are given in |E2j. 



IV. QUANTUM COMMUNICATION OVER LONG DISTANCES 
A. Concept of the quantum repeater: Nested entanglement purification 

We have now all necessary tools available to introduce the concept of the quantum repeater. Our goal is to create 
an EPR pair of high fidelity between two distant locations. Since nonlocal entanglement between distant partides 
cannot be created using only local operations, this involves the usage of a quantum channel, which is noisy in general. 
The bottleneck for communication over large distances is the sealing of the error probability with the length of the 
channel. Whcn using, for example, optical fibers and single photons as a quantum channel, both the absorption losses 
and the depolarization errors scale exponentially with the length of the channel. The state of the photon or the photon 
itself will therefore be destroyed with almost certainty if the channel is longer than a few half-lengths of the fiber. 

To overcome this limitation, we divide the long channel into N smaller segments and create less distant EPR pairs 
across each segment as visualized in Fig. The number of segments N is thereby chosen in such a way that it is 
possible to create EPR pairs with sufficicntly high initial fidelity F > F m j n ov er the distance of such a segment. In a 



next step, we connect these "elementary" pairs as described in Section III B. This leaves us with a pair between A 
and B with reduced fidelity Fn as given in (JÏ9|). In principle, one could now create many pairs between A and B in 
a similar way and then usc this ensemble of pairs for purification. But purification is only possible if the fidelity of 
the initial pairs is above a certain threshold value -Fmin as we have seen in Section III A[ This limits the number of 



of pairs one can connect before purification becomes impossible. We therefore connect a smaller number L -c N of 
pairs so that the resulting fidelity Fl stays above the threshold value for purification [Fl > F min ) and purification is 
possible. 

The general strategy will be to design an alternating sequence connection and (re-) purification procedures in such a 
way that the number of resources needed remains as small as possible, and in particular does not grow exponentially 
with TV and thus with l. In the remainder of this Section we describe a nested purification protocol which consists of 
connecting and purifying certain groups of pairs simultaneously in the following sense (see Fig. [ïï]). For simplicity, 
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assume that N = L n for some integer n. On the first level, we simultaneously connect the pairs (initial fidelity jFi) 
at all the checkpoints except at Cl, C2L, ■ ■ ■ , Cn-l- As a result, we have N/L pairs of length L (and fidelity Fl) 
between Ah Cl, Cl & C2L and so on. 

To purify these pairs, we need a certain number M of copies that we construct in parallcl fashion. [For keepin 
track of the resources, it is convenient to arrange them in form of an array of elementary pairs as is donc in Fig 
for L = 3 and M — 4]. We then use these copies on the segments A & Cl, Cl & C2L etc, to purify and obtain 
one pair of fidelity > F\ on each segment. This last condition determines the (average) number of copies M that we 
need, which will depend on the initial fidelity, the degradation of the fidelity under connections, and the efficiency 
of the purification protocol. The total number of elementary pairs we used up to this point is LM . On the second 
level, we connect L of these larger pairs at every checkpoint CkL (/c = 1,2...) except at Cl 2 , C2L 2 , ■ • ■ , Cn~l 2 - As a 
result, we have N/L 2 pairs of length L 2 between A & Cl 2 , Cl 2 & C 2 l 2 , an d so on of fidelity > Fl- Again, we need 
M parallel copies of these long pairs to re-purify up to a fidelity > F\ . The total number of elementary pairs involved 
up to this point is (LM) 2 . We iterate the procedure to higher and higher levels, until we reach the n-th level. As a 
result, we have obtained a final pair between A & B of length N and fidelity > F\. In this way, the total number R 
of elementary pairs will be (LM) n [where M n alone gives the number of required 'parallel channels' in Fig. |IÏ| . We 
can re-express this result in the form 

R = N l °^ M+1 (20) 

which shows that the resources grow polynomially with the distance N. 

The central feature of this nested purification protocol is the two-step process connection - purification on each 
nesting level. This purification loop is visualized in Fig. The curves shown in this figure are described by the 
formulas @ for purification (upper curve), and ( |Ï9|) for connection (lower curve, with TV replaced by L), respectively. 
Let us consider a given nesting level k, where we have N/L k ~ 1 pairs of fidelity F each. Starting from F, the fidelity 
Fl after connecting L pairs can be read off from the curve below the diagonal. Reflecting this value back to to the 
diagonal line, as indicated by the arrows in Fig. [ï^, sets the starting value for the purification curve. If Fl lies within 
the purification interval (-Fmm, F max ), then the purification curve lies above the diagonal, F' (Fl) > Fl, and iterative 
application of (^|) leads back to a fidelity larger or equal to the initial value F. Each step in the 'staircase' of Fig.|Ï2| 
corresponds to a successful purification step, and with help of the total number of necessary steps k max needed to 



recover one can calculate the needed resources M as described in Section III A . Once the initial value F of the fidelity 
is reobtained, we have N/L k pairs and we can start with the next level k+1. In summary, each level in the protocol 
corresponds to one cycle in Fig. ^2|. 

There are two conditions that have to be fulfilled to realize a closed loop. First, the fidelity after the connection has 
to be larger than the minimum value for purification (lower fixpoint), Fl > -F m in- If Fl < F m i n , then repurification 
is impossible. Second, the fidelity F one starts with and wants to reach again has to be smaller than the maximal 
reachable fidelity (upper fixpoint), F < F max . These two conditions determine both a threshold value for the error 
tolerances of the local operations and they limit the value L (the number of pairs which can be connected before 
purification becomes necessary). Please note that, being polynomial in L, the lower curve gets steeper and stceper 
near F = 1 for higher and higher vàlues of L. From this one can see that, for a given initial fidelity F, there is a 
maximum number of pairs one can connect before repurification becomes impossible. 

The threshold value for the nested purification scheme, i.e. the repeater, is tightcr than the threshold for simple 
purification on a single segment. Consider, as an example, the situation in Fig. [ï^. The chosen value p 2 = 0.97 of the 
error parameter (with rj = 1) lies above the purification threshold of 0.95 that we have found after Eq. (|Ï2|). On the 
other hand, for p 2 — 0.95, the purification interval in Fig. |Ï2] would shrink essentially to a point at F m j n = F max = 0.75. 
Clearly, in this situation no loop can be realized; even if we start initially at F = 0.75, the purification "interval" is left 
after the first connection process and the fidelity would subsequently converge towards the trivial fixpoint F = 0.25. 

When we employ the protocols A or B in the purification part of the nested scheme, we find error thresholds that 
are typically in the percent region. For scheme C, the situation is similar, as long as we do not use depolarization 
after each purification step (i.e. keep all involved pairs during the purification process in the Werner form). If we do 
use depolarization, thus realizing a modified version of protocol A, the loop cannot be closed. Differently speaking, 
for Werner states, the fidelity that is lost by the connection process can never be regained by purification with this 
scheme. For this to be true, the repurification condition F mSíX (FL(F)) > F must be fulfilled, where the upper fixpoint 
Fmax is a function of the fidelity of the auxiliary L-pair. It can be shown analytically that this is not possible for 
scheme C with Werner states. This underlines the fact that the purification scheme C is not a trivial variant of either 
of the schemes A and B. 

Note that both connection and purification are not smooth processes. This means that one can only perform whole 
steps and not parts of a step. Due to this fact, one will not get an exactly closed purification loop, in general, but the 
final fidelity after the purification (this is the value at which the ladder in Fig. O ends) may be slightly larger than 







the initial fidelity; secondly, the final fidelity may be different at each nesting level, but it is always larger or equal to 
the fidelity of the elementary pairs. 



B. Physical resources 

The length of a segment is limited by the transmission errors, and the number of segments N one splits up the total 
channel should be varied in order to find out the optimal configuration. The quantum repeater can be characterized 
by two quantities, which depend on the used connection and purification protocol 

• The total physical resources per segment needed to build up the EPR pair - this gives the number of necessary 
'paral·lel channels' between two checkpoints. 

• The total time which is needed to build up the EPR pair. 

The way to calculate these quantities is similar for schemes A and B, while the properties of scheme C are totally 
diffcrent. 



C. Physical resources per segment 

The physical resources per segment R se gment gives us the number of required connection lines between two check- 
points, and also half of the partides which should be stored at each checkpoint. On the borders (at A and B), it is 
equal to the number of partides to be stored, while at each checkpoint Ci , one has to take care about the incoming and 
outgoing string and therefore needs the double number of partides. Rsegment (vertical axis in Fig.^TJ) is completely 
determined by the resources M needed for each purification loop and is given by Rsegment = M n , where n = loq r,N is 
the number of nesting levels. The resources M needed for purification can be calculated as described in section III A . 

For schemes A and B, the physical resources are equal to the resources, since all pairs have to be created and 
therefore stored at the beginning of the process. Figure shows the resources M needed for a single purification 
loop, where it is assumed that L — 2, i.e. only two pairs are connected before re-purification takes place. M is plotted 
against a 'working fidelity', which is the fidelity one starts with before connection and one ends up after re-purification. 
Although the difference between schemes A and B does not look very dramàtic, note that M has to be t aken to the 



n power to calculate the resources per segment, where for long range communication n ss 10 (see Sec. IV E ). One 
can see that there exists an optimal working region, which can be undcrstood as follows: For small working fidelities 
F, two things happen: Firstly due to the polynomial law for the connection, the fidelity decreases strongly. Secondly, 
purification becomes less effective near the lower fixpoint Finin (the gain per purification step gets smaller). For both 
reasons, many purification steps are needed to recover and thus the needed resources are large. For working fidelities 
F close to the maximal reachable fidelity (upper fixpoint F max ), one does not loose much due to the connection, but 
purification is less effective - the gain per purification step is smaller. Therefore the the number of necessary steps is 
again quite large and so are the needed resources. In between there exists a region which is optimal in the sense that 
the needed resources are minimized. Here a working fidelity F sa 0.95 turns out to be optimal for the chosen error 
parameters. 

The error dependence of the resources is shown in Fig.[Ï4| . For smaller error probabilities (larger error parameters), 
two things happen: The optimal working fidelity gets larger and the needed resources M decrease. For imperfect 
local operations, the needed resources increase to infinity if one really wants to reach the maximal achievable fidelity 
-Fmaxj only for perfect operations and F = 1, this is not true because one starts with perfect pairs and does not loose 
due to connection. Therefore no purification is necessary and M = 1. 

Using scheme C for purification, the situation is different. The vertical axis of Figjïï] (which corresponds to the 
resources) is translated into a "temporal axis". Instead of creating all needed pairs at the beginning and operating 
parallelly on ensembles of pairs at each checkpoint (as in scheme A and B), one creates the pairs needed for purification 
every time one needs one and operates sequentially on the pairs. To perforin the purification process, every time only 
two pairs are involved. One is needed to store the purified state and the other one is repeatedly created to purify the 
first pair as visualized in Fig^. The number of physical resources required for purification is thus only 2, but every 
time the purification was not successful one has to start from the very beginning. Using scheme A or B one obtains 
a purified pair after k m ax successful purification steps, which is in this case also the number of really performed 
purification steps. Looking at Fig.^], one sees that a purification step which was not successful eliminates one branch 
of the tree, but after 4 steps (in this example) one ends up with a purified pair. In scheme C, on the other hand, one 
does not follow all branches parallel but follows the branches sequentially (and using always the same elementary pair 
for each purification step). This limits the number of involved partides, but increases the required time. 
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Looking now at the nested algorithm, it turns out that the required number of physical resources for scheme C grows 
linearly with the number of nesting levels and thus logarithmically with the number of segments (and the distance). 
This is because one additional pair is needed for storage purposes at each nesting level. One can understand this by 
inspecting Fig.^. First, three elementary pairs are connected (line 4) and used to purify the pair at line 3, which now 
is the 'elementary pair' for the next nesting level. Second, three of these re-purified pairs (line 3) are connected and 
used to purify the pair at line 2. Since one has to create such a pair repeatedly and all partides at lines 3 and 4 are 
involved, it is necessary to store the pair to be purified (here in line 2). Thcrcfore one needs one additional particle to 
store the pair at each new nesting level. This particle is not needed in all checkpoints, but only in those which lie on 
the borders of the corresponding nesting levels. That is why the maximum number of additional partides (physical 
resources) is required at the outermost places (A and B), while e.g. at checkpoint C±, no additional particle is needed. 
The maximal number of additional partides (physical resources) grows only linearly with the number of nesting levels 
and therefore logarithmically with the number of segments in contrast to the polynomially growing resources when 
using scheme A or B. The maximal resources per segment (in this case partides needed to store the pair) can be 
calculated by 

Rsegment = 1 + log L N = 1 + U (21) 



D. Temporal resources 

Here we will discuss the total time needed to create an EPR pair. This involves at least three parameters: 

• The typical time r op needed to perforin local operations (single qubit, two qubit and measurements). 

• The typical time T va j r to create an elementary pair. Using optical fibers and the model of the absorption free 
channel (AFC) pqffil , this can be expressed in terms of the other 2 parameters and is given by 

í '_3 Ss ment \ 

T P air = tafc = (5t op + 2T c i ass ) eV 'o / (22) 

where Iq is the half length of the fiber and hegment is the length of a segment. This just refleets the fact that for a 
single use of the AFC, all together 5 operations and 2 classical Communications (in this case the transmission of 
a photon) are necessary. The exponential function gives the average number of repetitions which are necessary 
due to absorption losses. 

• The time r c ; ass needed for classical communication to broadeast the result of the measurement. This time is 
determined by the length of the segment and the speed of light c and is given by 

^segment /'OQA 

Trfass — 

C 

To simplify the discussion, we will assume that the number of pairs L which are connected at each nesting level is 
some power of 2, L — 2 l , because in this case the connection process can be performed in parallel. For each connection 
process, 3 elementary operations (CNOT, measurement, operation depending on the outeome of the measurement) and 
classical communication over the distance of one segment are needed. For purification, also three elementary operations 
(e.g. for scheme A depolarization, bilateral CNOT and bilateral measurement) and classical communication of the 
result are necessary. 

Using scheme A or B, the time needed to perform a purification loop (the connection of L = 2 l pairs and re- 
purification, where k max successful purification steps are necessary) at nesting level m is given by 

tioo P (m) = 31t op + f(m)(2 k - l)r ci 

ass ~t~ fcmax 

(3T op + f{m)2 l T c i ass ) (24) 

where f(m) = (2 l ) m gives the length dependence of the elementary pair on the nesting level. The first two 
terms give the time required to perform the connection process in parallel fashion, where l connections and classical 
communication over the distance of (2 l — 1) times the length of an elementary pair is needed. The purification process 
involves k max purification steps and for each purification step classical communication over the distance of 2 l times 
the length of an elementary pair (since L = 2 l pairs werc connected before). One only has to take into account the 
successfull purification steps, since one runs through the whole purification tree (see Fig.||) in a parallel fashion. To 
calculate the total time needed to perform the nested algorithm, one can use the iteration formula 
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Uot{m) = t tot (m - 1) + ti oop {m) (25) 

where m corresponds to the m th nesting level and t to t (n) gives the total time needed for the creation of an EPR pair 
between A and B. To complete tiris iteration formula, note that t to t(0) = T pa i r which is just the time needed to create 
an elementary pair at the beginning. The total time if one has N — L n segments and L = 2 l is thus 

^tot ~ ^(3^ ~1~ 3fc maa; )T p + ^2 1 + k max 2 í - Ttfass -\- Tpair (^) 

The main contribution is given by the classical communication if r op is not too large and is just t tot w (2k max + 

l)NT class . 

Using scheme C for purification, the time needed to perform a purification loop on nesting level m is given by 

tioop{m) = (max{í ioop (m - 1)} + iW op + f{m)(2 l - l)r c/ass ) M + S(3r op + f (m)2 l T class ) (27) 

where the first term is the time needed to create a pair for purification (by building up 2 l pairs at nesting level 
(m-1) and connecting them), which must be done sequentially M times. The second term gives the time needed for 
purification (operations and classical communication over the distance of 2 l times the length of an elementary pair) , 
where S purification steps have to be performed. The function f(m) — (2') reflects the length dependence of the 
elementary pair on the nesting level and the time needed to create an elementary pair at the lowest nesting level 
T pa ir = í/oop(0), which completes this iteration formula. The maximum appearing in this expression is taken over 
all 2 l segments on nesting level (m — 1). This is necessary because all 2 l new elementary pairs have to be built up 
before the next connection process can start. Therefore one has to wait for the slowest one. Note that the time for 
purification is fixed for schemes A and B (obtaining bad results in this case corresponds to the elimination of a branch 
in Fig.|^) , while for scheme C the purification process has to be started from the very beginning when obtaining a bad 
result. This results in the fact that the performed purification steps S and the used purification resources M (and 
therefore the needed time) might vary from case to case. This formula is sufircient to calculate the total time needed 
to perform the nested algorithm up to nesting level m and t tot (n) = ti oop (n) gives the total time needed to build up 
an EPR pair between A and B. 

A rough estimation of the total time can be obtained when using (15) and ( fÏ7| ) to calculate the average number of 



purification resources M and purification steps S. In this case, one can write down a closed expression for the total 
time, which is polynomially growing with the number of segments N and thus with the distance. Comparing this 
time with the results obtained from a simulation of the whole process, it turns out to be smaller by a factor of ~ 3. 
In some sense, the vertical axis (physical resources) in Fig.|TT| is translated into a temporal axis. 



E. Comparison of the 3 schemes 



We will give some typical numbers to quantify the properties of the nested purification protocol if the purification 
loop is performed with hclp of scheme A,B or C. The following tables are based on error parameters of = p% = 

i] = 0.995) and a working fidelity F=0.96, the fidelity one starts with at the lowest nesting level and one ends up 
after each purification loop. The length of a segment is assumed to be in the order of half length of an optical fiber, 
Isegment ~ Wkm. In this case , the number of segments is proportional to the distance between two locations. If the 
number of segments is in the order of N = 2 7 = 128, this will be called 'Continental scale' (distance of m lOOOfcm), 
while N = 2 10 = 1024 is called 'Intercontinental scale' (distance of rí lOOOOfcm), since one can create an EPR pair 
between two international cities. It is assumed that the connection is performed in paral·lel. The columns from left 
to right have the following meaning: 

• resources: this gives the physical resources per segment. 

• time: gives the total time in seconds needed to create an EPR pair. The number is based on the application 
of the AFC to create elementary pairs (r pa ir — tafc — 3.2 * 10~ 4 s). The half length of the fiber is considered 
to be lOkm, which is also the length of the elementary segment (r c / ass = 0.33 * 10 _4 s). All operations (single 
qubit, 2 qubit, measurement) are considered to be performed in r op = 10~ 5 s. The time is calculated using a 
simulation of the nested entanglement purification process, averaging over a few hundred runs. 





Continental scale 


Intercontinental scale 




resources 


time 


resources 


time 


A 


1.58* 10 y 


3.88* lO" 2 


9.01 * 10 12 


0.298 


B 


329 


1.34* 10"^ 


4118 


0.103 


C 


7 


0.77 


10 


15.69 
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One sees that the physical resources needed when using scheme C are some orders of magnitude smaller, while the 
time needed is larger. For a practical implementation, the achievable bit rates are not very impressive, but note that 
all numbers get much better if the local operations can be performed better (error parameters closer to 1) as can be 
seen in FigJÏJ. 

So the reachable fidelity is larger and at the same time, the required resources are smaller when using better local 
operations. Note that the total time needed is in the order of the time for classical communication over this distance, 
which is just t class = 3 ffgt%w„ = °- 0Ms - 



V. SUMMARY 



We have shown that it is possible to create EPR pairs via a noisy channel, with an overhead in physical resources that 
grows only logarithmically with the lcngth of the channel, while the time needed for the crcation grows polynomially. 
The central idea is to use a nested purification protocol for connecting a sequence of EPR pairs at certain "connection 
points" within the channel, whose roles are reminiscent of to classical repeaters. Different from the classical situation, 
the concept of the quantum repeater is not a local amplifier, but it involves both the local checkpoints and global 
purification protocol. Our scheme tolerates errors for local operations and measurements that are in the percent 
region. 

This work was supported in part by the Osterreichischer Fonds zur Fòrderung der wissenschaftlichen Forschung, 
and by the European TMR network ERB-FMRX-CT96-0087. 



Permanent address: Institut für Theoretische Physik, Universitat München, Theresienstrasse 37, D-80333 München. 
B. Schumacher, Phys. Rev. A, 54, 2614 (1996); B. Schumacher and M. D. Westmoreland, ibid. 56, 131 (1997). 
W. K. Wootters and W. H. Zurek, Nature 299, 802 (1982). 
D. Dieks, Phys. Lett. A 92 271 (1982). 

R. J. Glauber, in Frontïers in Quantum Òptics, (eds. E. R. Pike & S. Sarkar), pp. 534-58 2. Bristol: Adam H l·lger. 
W. Tittel, J. Brendel, B. Gisin, T. Herzog, H. Zbinden, and N. Gisin, Los Alamos report ^uant-ph/9707042. 
For opti cal fibers, this leng th is typically lOkm; see e.g. ||. 

A. M. Steane Phys. Rev. Lett. 78 2252 (1997). D. Gottesman, |quant-ph/970202S| ; For a review 



P. Shor, |quant-ph/9605011 
see e.g. J. Preskill, |quant-ph/971204£ 



E. Knill and R. Laflamme, quant/ph-9608012. See also D. Aharonov and M. Ben-Or, ^uant-ph/9611025; C. Zalka, quant 



ph/ 961202S 



Rev. Lett. 77 , 2818 (1996). 

Wootters, Phys. Rev. Lett. 70, 1895 (1993). 



C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, J. A. Smolin, and W. K. Wootters, Phys. Rev. Lett. 76, 722 
(1996). 

C. H. Bennett, D. V. DiVincenzo, J. A. Smolin, and W. K. Wootters Phys. Rev. A 54, 3824 (1996). 
N. Gisin, Phys. Lett. A 210 151 (1996). 

D. Deutsch, A. Ekert, C. Macchiavello, S. Popescu, and A. Sanpera, Phys. 
C. H. Bennett, C. H. Brassard, C. Crepeau, R. Josza, A. Peres, and W. K. 
A. K. Ekert, Phys. Rev. Lett. 70, 661-663 (1991). 

G. Giedke et al., unpublished. 

H. J. Briegel, W. Dür, J. I. Cirac, and P. Zoller, quant-ph-9803056 (1998). 
C. Miquel, J. P. Paz and W.H. Zurek, Phys. Rev. Lett. 78, 3971 (1997). 

A similar form as in (Q) may be obtained for two-bit rotations with a corresponding isotropic noise function, although the 
calculation is more complicated. In particular, we have not yet been successful in giving a proof in terms of a manifest 
parametrization-independent average over the invariant group manifold of SU(4). Isotropic noise in this case means that 
the depolarizing channel does not distinguish between the two qubits and acts on them in a symmetric way. 
[19] This assumption is not essential; it is motivated by the fact that the depolarization is a itself an inessential although 
simplifying step in the Bennett et al. protocol. If the reliability p\ is the same for all one-bit operations involved in the 
realization of the twirl, an imperfect depolarization can effectively be described by 



Tp 12 = PTT ld p' 12 + 



Pt 



Il2 



where pr = (pi) 2 summarizes reliability of the total twirl operation T and p'12 is the bilocal state before the twirl is applied. 



13 



This operation still depolarizes p' 12 and brings it to the Werner form. However, the fidelity F" of the resulting state is 
reduced with respect to the fidelity F' of the state p' 12 in ([]) by the relation 

f » =PtF > + L=J!L. 

4 

[20] C. W. Helstrom, Quantum detection and estimation theory, Associate Press, London 1976. 

[21] Any sequence of operations at the same node, e.g. two single qubit operations Oi <g> O2 followed by a CNOT operation will 
be described by a single two-qubit oper ation O12 and consequently with only a single err or parameter 



[22] W. Dür, Master Thesis, available under http://www.uibk.ac.at/c/c7/c705/staff/wd.html 
[23] see also N. Lütkenhaus, unpublished. 

[24] A. Zeilinger et al., Phys. Rev. Lett. 78, 3031 (1997) In their scheme, however, a three-particle entangled state is obtained 

rather than a distant EPR pair. 
[25] This procedure is not necessary, 

all operations in terms of (pH7|). 
[26] S. J. van Enk, J. I. Cirac, and P. Zoller, Science, 279, 205 (1998). 
[27] S. J. van Enk, J. I. Cirac, and P. Zoller, Phys. Rev. Lett, 78, 4293 (1997). 



1 2 

/ O-a rO \ 

/ I I 1 

11 II 

vO-7 ^> , 



4 

FIG. 1. Partides required for purification in Standard recurrence schemes. 
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FIG. 2. Typical purification process. At each purification step, either both pairs are discarded (if the purification was not 
successful) or one pair is discarded (if the purification was successful). The left-over pairs are again used for purification at the 
next step. 
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FIG. 4. fidelity F plotted against the number of successful purification steps for scheme B (upper curve) and scheme A. 
Based on fixed error parameters P2 — rj of 1% and initial fidelity Fo = 0.7. 
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of scheme A (x) and scheme B (•) plotted against error parameters P2 = rj. 



FIG. 6. Typical purification process. At each step, the same elementary pair no is used. If one purification step is not 
successful, one has to start from the beginning. 
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FIG. 7. Reachable fixpoint plotted against parameterization parameter e (shape of the state) for fixed initial fidelity Fo = 0.7. 
Plots from bottom to top correspond to error parameters p2 = rj of 4%, 3%, 2%, 1%, 0%. 
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FIG. 8. Reachable fixpoint plotted against initial fidelity Fo for scheme B (constant function) and scheme C. Based on fixcd 
error parameters p 2 = r\ of 4% and binary pairs (A,C). 
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FIG. 9. Partides involved in the connection process. Particle 1 is located at A, partides 2 and 3 at Ci and partide 4 at B. 
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FIG. 10. Connection of a sequence of TV EPR pairs. 




FIG. 11. Nested purification with an array of elementary EPR pairs. 



17 



1 



± max 


















- 
















/ F L 




0.3 0.4 


0.5 


0.6 


0.7 0.8 


0.9 


1 



F 

FIG. 12. "Purification loop" for connecting and purifying EPR pairs. The parameters are L = 3, 77 = pi = 1, and P2 = 0.97. 
The (upper) purification curve corresponds to scheme A, Eq. (S). The (lower) connection curve is described by (B with N = 3. 
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FIG. 13. Resources M for purification versus working fidelity F for scheme A (upper curve) and the scheme B (lower curve). 
The errors of all operations are 0.5%, and L — 2. 
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FIG. 15. The nested purification algorithm using the Innsbruck protocol. At each nesting level, one additional particle at 
the borders of tbis nesting level is needed to store the pair meanwhile 
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FIG. 16. Additional partides needed at each checkpoint for the Innsbruck protocol. The number of segments N — 2 4 = 16, 
L—2 pairs are connected at each nesting level and the number of nesting levels is n=4. 
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